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Abstract: We present four results for oscillons in classical (j)'^ theory in D+l space-time 
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oscillons before, during and after its ultimate demise. 
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1 Introduction 



Oscillons are non-top ological, long-lived and quasi-periodic field configurations in certain 
non- linear field theories [1-12]. Much is known both analytically and numerically about 
these objects, but a complete understanding of their existence, creation and ultimate decay 
is still not settled. 

In the following, we will give at least partial answers to four questions one may want 
to ask about oscillons in general, but concentrating on a real, classical scalar (/)^-theory. 
The action reads 



S 



A 



(1.1) 



and it is known that the theory has oscillon solutions. The two degenerate minima are at 
V = ib/i/\/A and the mass of excitations around these minima (the "radiation" frequency) 
IS m = V2fi. We will assume spherical symmetry in D dimensions, and in this case the 
classical equation of motion reads: 



D - I 



dr 



0. 



fl.2) 



We will discretize and solve this equation on a spatial lattice in real-time. By assuming 
spherical symmetry, the lattice is in practice one-dimensional, with the physical dimen- 
sionality appearing only as the parameter D. This being the case, we will treat D as a 
continuous parameter, i.e. allow for a non-integer number of dimensions. 

In order that radiation emitted from the oscillon does not influence the later evolution 
of the oscillon, we will impose absorbing boundary conditions, generalizing the approach 
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in [3] to D dimensions. This means that at the large-r boundary, we replace (1.2) by (see 
Appendix A) 



+ drdt + ^ ^ ^ dt + /i^ 
Zr 



0. (1.3) 



The alternatives are to have periodic boundaries and large lattices or some dissipative 
procedure far from the oscillon core [4, 5]. At r = 0, spherical symmetry imposes drff) = 0, 
and we therefore have 

- Ddl - ^2 ^ A (A^] (/> = 0. (1.4) 
We initialize the field with a Gaussian profile 

(/.(r,0) = 1-Cexp(^-^^ , (1.5) 

with two free parameters, the amplitude C and the width tq. The initial momentum is set 
to (j){r^ t = 0) = 0, for all r. 

We will take n = 1 and A = 1, which amounts to a rescaling of field and space-time 
coordinates, so that t— ?• fit, x— ?• ^x, (j) ^ (j)/v. Then the minimum of the potential is at 
= ±1, and lifetimes and sizes of the oscillons are in mass units throughout. The radiation 
frequency is then m = \/2. 

In the following, we will consider the evolution of the energy, as well as the instanta- 
neous oscillation frequency of the center of the oscillon uj{t), taken to be 

u{t) = — , (1.6) 

where 6t is the time between a crossing of 0(0, t) = and the crossing a period (i.e. two 
crossings) later. 
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Figure 2. The lifetime of oscillons from Gaussian initial conditions in = 3, in the C-rg plane. 



Fig. 1 (left) is an example of the energy inside a box of radius r = 15 around the 
oscillon (red), the radiation in the range 20 < r < 90 (blue) and the sum over the whole 
lattice (magenta). At first the oscillon sheds energy into radiation (red curve drops, blue 
curve rises), but total energy is conserved. Then around t = 100, the radiation reaches the 
absorbing boundary, and radiation and total energy drops, while the oscillon continues to 
pump energy into radiation. This is in fact only a transient stage until the Gaussian initial 
condition settles to a true oscillon at late times, with some almost constant energy. 

Fig. 1 (right) is the later time evolution of the oscillon, we see the brief transient and 
the long oscillon stage, which ends suddenly at i ~ 7200, where all the energy is dumped 
into radiation. We also note that during the oscillon stage, very little radiation is emitted, 
ad so very little radiation reaches the boundary. 

2 What is the hfetime of the oscillon in D dimensions? 

In [4], spherically symmetric simulations in D+1 dimensions were performed using dis- 
sipative dynamics away from the oscillon. Using absorbing boundary conditions allows 
much more efficient numerics, since we can do with smaller lattices. Also we find that this 
approach is better at getting rid of spurious waves bouncing back from the far boundary. 

In Fig. 2 we show the lifetime of the oscillon in D = 3 as a function of the width and 
amplitude of the initial Gaussian. The lifetime is simply measured from the initial profile 
at t = until the oscillon collapses. We see two clear regions where the lifetime is large. 
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Figure 3. The lifetime of oscillons from Gaussian initial conditions in Z? = 4 (left) and D = 5 
(right), in the C-rg plane. 



up to t = 8000, corresponding to starting the oscillon either up the potential (left side) or 
across to the other potential minimum and beyond (right side). The two parameter regions 
are of course connected by the oscillating behaviour of the oscillon, so that they should not 
really be seen as separate. 

The black region is where the initial condition does not evolve into an oscillon, and 
it is clear that there is a lower limit in the amplitude, around the inflection point of the 
potential C > 1 — \/l/3 ~ 0.42. Although more fuzzy, there also seems to be upper and 
lower limits to the region of large lifetime, say around —2 and 4. The limits on the allowed 
width is strongly correlated with the amplitude, but in general ro > 2. Finally, there may 
be an upper limit to the width of roughly tq < 5, but this is less clear. 

In Fig. 3 we show a similar result, but now for D = A and D = 5. We observe a 
similar picture of two main parameter regions, but whereas the lower limits to the basin of 
attraction are again quite sharp, the upper limits are more fuzzy. As in D = 3, the lower 
limits are ro > 2, and the amplitude above the inflection point, larger for small width. But 
the maximum lifetime regions are at larger width and amplitude than for D = 3. 

Perhaps surprisingly, the maximum lifetime in D = 5 is larger than for D = 4. In 
fact, because the lower limits shift, and the maximum lifetime regions move about, a given 
parameter set C, ro can fall out of the oscillon basin of attraction and into the collapse 
region (black area). In [4], the lifetime as a function of D was reported to follow a power 
law until D = 5, after which a second power law set in. In our language, the configuration 
considered there is (C, ro) = (2,3.3), which for D = 3 is nicely in the maximum lifetime 
region, but which for D = 4, 5 falls off the basin of attraction. In [4] the onset of the second 
power law was used to argue that oscillons do not exist for high D. More on this in section 
4. 

In Fig. 4 we show the C-ro plane for D = 2. It is known that oscillons can live 
for millions of oscillations (see for instance [6]), and we simply cut off at t = 10000, 
since we are interested in the structure of the phase space. We see that indeed two major 
regions establish themselve, with very long lifetimes, but now on the left-hand side, multiple 
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Figure 4. The lifetime of oscillons from Gaussian initial conditions in = 2, in the C-rg plane. 

"ridges" areopening up. Also for large enough amplitude, oscillons arise from very narrow 
initial Gaussians (small tq). Again, from the point of view of a fixed parameter point (say 
(C, ro) = (2,3.3)), the right-hand oscillon region has shifted down and left. 

As was beautfifully demonstrated in [8, 9], the lifetime as a function of width and 
amplitude in fact has a spiky structure, where by carefully tuning the parameters, one can 
get perhaps arbitrarily large lifetimes. In Fig. 5, we show with much higher resolution 
the region near the point tq = 2.28, C = 2 in D = 3, which is one instance of the spiky 
structure reported in [8, 9]. In fact, the substructure is a set of very narrow bands, and 
indeed factors of 3 or more in lifetime can be achieved by careful tuning. As we see, the 
tuning has to be extremely fine (here less than one per mille) to see the additional structure. 
The broader behaviour of the lifetime is of course represented by the coarser sampling of 
figures 2-4, and a population of oscillons created during a phase transition are likely to 
have these lifetimes [13-16] , but the fine structure is a signal of the non-trivial nature of 
the dynamics, and deserves closer scrutiny in its own right. 

3 Is there more than one oscillon? 

The current understanding of oscillons is that they are attractor solutions to the equations 
of motion, in the sense that within a basin of attraction, all initial localized field config- 
urations approach some oscillon solution, we have seen in the previous section how these 
regions manifest themselves in terms of the lifetime. 
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Figure 5. The lifetime in D — 3 in a narrow, but highly resolved region the C-tq plane, showing 
a fine band structure. 
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Figure 6. A number of trajectories of uj{t) starting from different Gaussian initial conditions in 
-D = 3. Left; Initial conditions leading to oscillons. Right: Initial conditions in the black region, 
where there is immediate collapse to radiation. 



The approach to a true oscillon happens by an initial shedding of energy as radiation 
and the changing of shape to the oscillon profile. The oscillon performs a non-linear 
oscillation with a slowly changing frequency uj{t) slightly below the mass m = V^/j., which 
is why it cannot easily decay into free excitations of the field. 

Oscillons slowly evolve by emitting a little radiation and changing their frequency to 
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a critical value Wc, where they collapses into only radiation at uj = m (see for instance [4]). 
It is therefore natural to think of an oscillon as a whole trajectory of field configurations 
ending at ojc, and different initial configurations in the basin of attraction will asymptote 
to different points along this trajectory. 

Now the question is: Is there more than one oscillon trajectory in (j)'^ theory in D + 1 
dimension? In other words, can the field space (of localized energy configurations) be 
split up into only two basins of attraction, one leading to the oscillon (trajectory) and one 
leading to immediate collapse into radiation? 

The first hint that the oscillon is unique is that the basins of attraction are dominated 
by two large regions, corresponding to the extremes of the oscillation (Figs. 2-4). Fig. 6 
(left) shows the trajectories of uj{t) of the oscillon center (j){0,t), for a number of Gaussian 
initial conditions, shifted so that the final collapse is simultaneous. The initial conditions 
are picked from different mutually disconnected regions in D = 3. We see how they 
asymptote either to the same, unique trajectory, for the short-lived with an additional 
oscillation due to the oscillon radius "wobbling" . In particular the evolution after collapse 
is remarkably similar. 

In Fig. 6 (right) we picked a number of points from the black collapse region, and we 
see a completely different behaviour. u}{t) > m throughout, and it is not even clear that 
there is a well-defined oscillation. At late times the co again asymptotes to m. Obviously, 
we cannot claim that gaussian initial conditions sample all configuration space, but the 
results suggest that the oscillon trajectory is unique, and that in this sense, there is only 
one oscillon. We checked that a similar picture emerges in L) = 2,4 and 5. 

4 Can there be oscillons in a large number of dimensions? 

In [7], it is suggested that oscillons cannot exist for a large number of spatial dimensions. 
This is based on analytical calculations for gaussian profiles, which show that there is a 
minimum Gaussian width ac{D) which allows for oscillons. 



and since this diverges around D ~ 6.88, no oscillons can exist for D = 7 and beyond. 

As shown in Fig. 7, at least the value of this limit must be put down to restricting 
not only the initial condition, but the whole evolution to a Gaussian anzats. We show the 
minimum Gaussian width leading to an oscillon at C = 2, i.e. the smallest width along 
a particular slice in C-ro space, as a function of D. We see that although this quantity 
increases with D, there is no sign of any divergence near D = 7. We have superposed 
the Gaussian analytic result of [7] and a linear fit. It is however true that because the 
lifetime becomes smaller and smaller, ultimately it is hard to decide whether the oscillon 
configuration has the time to establish itself, before it collapses. But this is an artifact of 
the Gaussian initial condition, rather than the absence of a quasi-breather solution (see for 
instance [4]). 
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Figure 7. The minimum width r„iin allowing for an oscillon at fixed amplitude C = 2. Superposed, 
the Gaussian estimate (magenta) and a linear fit rmin = (0.512 ± 0.004)D + (0.71 ± 0.01) (blue). 




Figure 8. The lifetime of oscillons from Gaussian initial conditions in I? = 6 (left) and D = 7 
(right), in the C-rg plane. 

Fig. 8 shows the lifetime in the C-rQ plane for D = 6 (left) and D = 7 (right). We see 
that there are certainly oscillons of significant lifetime. We should note that the left-out 
regions in the upper left and right-hand corners suffer from numerical instabilities which 
we were unable to cure. This becoms worse at higher D, and prevents us from simulating 
beyond D = 7. But the two bulk regions and their lower boundary to the no-oscillon basin 
are physical. In D = 7, these are rather narrow, and it is possible that there is an upper 
limit to D, or that only highly tuned bands remain. 
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5 What happens when an oscillon dies? 

At the end of its life, the oscillon suddenly collapses and all its energy is emitted into 
radiation, free field fluctuations peaked at uj = m. As was also remarked in [3], the non- 
linear oscillation of the oscillon can be resolved by Fourier transform into a discrete set 
of frequency peaks at integer multiples of the basic frequency. And so although the basic 
frequency is below the mass threshold, the higher modes are not. These modes are in fact 
responsible for the slow emission of radiation during the oscillon stage, as was shown in 
[3, 4] (see also [12]). But because the secondary peaks are exponentially suppressed, very 
little radiation is emitted. 

In Fig. 9 (left) we first show the Fourier transform in time of the center of the oscillon, 
r = 0, in a time interval around t = 4000. The initial condition is the same as for 
Fig. 1. We see that the dominant frequency is indeed slightly below the radiation frequency, 
ojo/m = 0.955, and that there are exponentially suppressed peaks at integer multiples of 
this. Similarly, in Fig. 9 (right), we show the frequencies present far into the radiation 
domain, r = 80. Now the largest peak belongs to the 2 x wo-mode, whereas the basic 
frequency is present only through the part of its peak that reaches beyond uj/m = 1. And 
so the width of the peaks is important, corresponding to the oscillation not being exactly 
periodic. 

In Fig. 10 (left), we show the spatial (r-space) Fourier transform in the radiation 
region 20 < r < 90. The peaks in k/m correspond exactly to the uj in Fig. 10 through 
= m? + k'^, and we again observe the dominant 2a;o-peak and the smaller cjo-tail near 
k = 0. We can now track the amount of radiation in each of these peaks (near 2uJo and m) 
by integrating up the area under each peak and displaying their evolution in time. This is 
shown in Fig. 10 (right). We see that deep in the oscillon regime, the 2ujq peak dominates 
the radiation. But as early as t = 7150, the radiation peak takes over. This precedes the 
collapse of the oscillon by ~ 100 in time units, as can be seen from the energy in Fig. 1. 

This means that the death of the oscillon is not an instantaneous process, but is 
preceded by a change in the radiation pattern. Oscillons decay because the width of the 



-9- 




primary frequency peak near uq eventually opens up enough that energy can be emitted 
through it. 



6 Conclusion 

We have considered four separate, but related, features of the oscillon in classical i;^^-theory 
in D dimensions. We have seen that 1) the oscillon is likely to be unique in the sense of there 
being only a single oscillon trajectory in w-space, or one non-trivial basin of attraction; 
2) we have computed the lifetime for Gaussian initial conditions in D = 2 — 7 and shows 
a two-region basin-of-attraction structure, the boundaries and maxima of which however 
shift as a function of D. Highly tunable lifetime "spikes" are in fact bands, in extremely 
small regions of parameter space; 3) there is no sign that oscillons cannot exist in D > 7, 
or that their lifetime should be short. In fact the maximum lifetime as a function of D 
has a minimum around D = 4; 4) the eventual collapse of the oscillon is preceded and 
triggered by the primary oscillation frequency becoming efficient at radiating away energy. 
And therefore not only the position of the peak in the spectral function but also the width 
is of importance for the stability of oscillons. 
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A Absorbing boundary conditions in D dimensions 

We consider the field far from the center of the oscillon, in terms of the displacement away 
from the minimum, ip = (j) — 1. The linearized equation of motion reads 



d-d' 



D - 1 



dr + m 



99 = 0, 



m 



(A.l) 
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Separating variables using if{r,t) = R{r) exp{—iojt) allows us to write (A.l) 

[r^d^ + {D- l)rdr + (w^ - m^y] = 0. (A.2) 
Using further R{r) = A'^~^^^'^S{r) and x = kr = ^ uP- — m^, we get the Bessel equation 

x'^S" ^S' [x^ -hj)-2f^S = ^. (A.3) 

In the limit r ^ 1, x ^ 1 this has the solution x"^/^ exp(ix), so that finally 

<^(r, t) = r(^-^)/2 [^(^^ _ ^ 

We now note that the outgoing wave w > obeys 

(^dr -ik + ^^^^ <f = 0. (A.5) 

Multiplying this by u), expanding k to first order in rv? juP' and making the replacement 
—iu: —7- c?j, we arrive at (1.3). 

References 

[1] M. Glciser, Phys. Rev. D 49 (1994) 2978 [hep-ph/9308279]. 

[2] E. J. Copcland, M. Glciser and H. -R. MuUcr, Phys. Rev. D 52 (1995) 1920 
[hep-ph/9503217]. 

[3] P. Salmi and M. Hindmarsh, Phys. Rev. D 85 (2012) 085033 [arXiv:1201.1934 [hep-th]]. 

[4] P. M. Saffin and A. Tranberg, JHEP 0701 (2007) 030 [hep-th/0610191]. 

[5] M. Gleiser and A. Sornborger, Phys. Rev. E 62 (2000) 1368 [patt-sol/9909002]. 

[6] M. Hindmarsh and P. Salmi, Phys. Rev. D 74 (2006) 105005 [hep-th/0606016]. 

[7] M. Gleiser, Phys. Lett. B 600 (2004) 126 [hep-th/0408221]. 

[8] E. P. Honda and M. W. Choptuik, Phys. Rev. D 65 (2002) 084037 [hep-ph/0110065]. 

[9] G. Fodor, P. Forgacs, P. Grandclement and I. Raez, Phys. Rev. D 74 (2006) 124003 
[hep-th/0609023]. 

[10] M. Gleiser and D. Sicilia, Phys. Rev. D 80 (2009) 125037 [arXiv:0910.5922 [hep-th]]. 

[11] M. Gleiser and D. Sicilia, Phys. Rev. Lett. 101 (2008) 011602 [arXiv:0804.0791 [hep-th]]. 

[12] G. Fodor, P. Forgacs, Z. Horvath and M. Mezei, Phys. Rev. D 79 (2009) 065002 
[arXiv:0812.1919 [hep-th]]. 

[13] M. Broadhead and J. McDonald, Phys. Rev. D 72 (2005) 043519 [hep-ph/0503081]. 

[14] M. Gleiser, N. Graham and N. Stamatopoulos, Phys. Rev. D 83 (2011) 096010 
[arXiv:1103.1911 [hep-th]]. 

[15] M. A. Amin, R. Easther and H. Finkel, JCAP 1012 (2010) 001 [arXiv:1009.2505 
[astro-ph.CO]]. 

[16] M. A. Amin, R. Easther, H. Finkel, R. Flauger and M. P. Hertzberg, Phys. Rev. Lett. 108 
(2012) 241302 [arXiv: 1106.3335 [astro-ph.CO]]. 



- 11 - 



